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ABSTRACT 

Dust emission at sub-millimetre wavelengths allows us to trace the early phases of star 
formation in the Universe. In order to understand the physical processes involved in 
this mode of star formation, it is essential to gain knowledge about the dark matter 
structures - most importantly their masses - that sub-millimetre galaxies live in. Here 
we use the magnification effect of gravitational lensing to determine the average mass 
and dust content of sub-millimetre galaxies with 250 /im flux densities of 6*250 > 15 mJy 
selected using data from the Herschel Multi-tiered Extragalactic Survey. The positions 
of hundreds of sub-millimetre foreground lenses are cross-correlated with the positions 
of background Lyman-break galaxies at 2; ~ 3 — 5 selected using optical data from 
the Canada-France Hawaii Telescope Legacy Survey. We detect a cross-correlation 
signal at the 7-a level over a sky area of one square degree, with ^ 80% of this signal 
being due to magnification, whereas the remaining ^ 20% comes from dust extinction. 
Adopting some simple assumptions for the dark matter and dust profiles and the 
redshift distribution enables us to estimate the average mass of the halos hosting the 
sub-millimetre galaxies to be log;^Q [Af2oo/-^0] = 13-17to o8(s^^*'-) ^^'^ their average 
dust mass fraction (at radii of > lOkpc) to be Mdust/-M200 ~ 6 x 10~^. This supports 
the picture that sub-millimetre galaxies are dusty, forming stars at a high rate, reside 
in massive group-sized halos, and are a crucial phase in the assembly and evolution of 
structure in the Universe. 

Key words: Submillimeter: galaxies. Gravitational lensing: weak. Galaxies: high- 
redshift 



1 INTRODUCTION 

In the current picture galaxies form in the centres of the 
gravitational potential wells of dark matter halos, when gas 
cools and star formation sets in ( [Eggenet al.|1962| ). However, 
many of the observed properties of galaxies - e.g., scaling 
laws between different observables like the IFaber fc Jack-I 
[son] ( |1976| ) relatio n the|Tully fc Fisher] (1 1977| relation, and 
the m-(j relation (Ferrarese & Merritt 20001) - remain to 



be explained by a complete theory of galaxy formation and 
evolution. We know from observations in the z < 1 Universe 
that star formation efficiency, i.e., the fraction of baryonic 
meiss turned into stars, depends on the environment, most 
importantly on the mass of the dark matter halo as shown 



by Leauthaud et al. ( 2012 1. The bulk of star formation, how- 



ever, happened at earlier cosmic epochs (Madau et al. 19961, 



when the Universe was less than half of its current age, which 
necessitates studies at high redshift (2 > 1). 

Young stars are still enveloped in clouds of gas and dust, 
and therefore most of their radiation is absorbed by dust, 
and re-emitted at far-infrared wavelengths. Since this is the 
dominant source of emission at these wavelengths, the lu- 
minosity of a galaxy in the far-infrared wavelength range is 
directly related to the amount of young stars, and therefore 
the rate of star formation. Sub-millimetre telescopes are sen- 
sitive to this radiation. Moreover, the detection of this light 
is less sensitive to redshift than in the optical, due to the for- 
tunate shape of the spectral energy distributions of galaxies 
in the sub-millimetre regime yielding samples of star- forming 
galaxies over a wide redshift range. Establishing a picture of 
the dark matter environment of these sub-millimetre galax- 
ies through measuring their total (i.e., baryonic plus dark 
matter) masses is hence a fundamental ingredient for under- 
standing the build-up of stellar mass over cosmic time. 

Gravitational lensing is the most direct method of mea- 
suring mass in the distant Universe, irrespective of whether 



it consists of dark or baryonic matter, and independent of 
any assumptions about its dynamical state. Thus, it repre- 
sents a powerful tool for measuring the masses of dark mat- 
ter structures, both using fewer astrophysical assumptions 
and being complementary to, e.g., velocity dispersion mea- 
surements. In many cases the lensing effect of an individual 
deflector is too weak to be detected. In weak gravitational 
lensing (WL; see e.g. |Bartelmann fc Schneid er|200l| a statis- 
tical approach is used, averaging the signals of many lenses 
and/or sources. 



In order to detect the effects of WL by a lens popula- 
tion, a suitable source population is needed that lies behind 
the lenses (from the observer's point of view). The extended 
redshift distribution of sub-millimetre galaxies means that 
most potential background galaxies are spatially unresolved 
in even the best ground-based optical data. Hence the tra- 
ditional shear technique of WL, which requires ellipticity 
estimates, cannot be used in this particular case, and one 
has to turn to the magnihcation technique, which does not 
require resolved sources ( [Van Waerbeke et al.|20Tol ). 



Here we measure the WL-magnification effect of a sub- 
millimetre galaxy sample to estimate their average dark 
matter halo density profiles. As background sources we 
use optically-selected Lyman-break galaxies (LBGs; |Steidel| 
|et al. l 1996) which are located at even higher redshifts. The 
sub-millimetre as well as the optical data are described in 
Sect. [2] An outline of the magnification technique is given 
in Sect. [3] Results are presented in Sect. [5] and discussed in 
Sect.|6] A summary and an outlook are give n in Sect.[7| We 
assume a WMAP7 cosmology throughout ( [Komatsu et al.| 
|20TT|): (^M, ^A, h, as) = (0.27, 0.73, 0.70, 0.81). 
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2 DATA 



2.1 Sub-millimetre data 



The sub-millimetre galaxies are selected from the Herschel 
Multi-tiered Extragalactic Survey (HerMES; [Oliver et al. | 
|2010[|HerMES Collaboration et al.|2012[ ), observed wit h the 
Spectral and Photometric Imaging Receiver (SPIRE; Grif- 
|fin et ar]|2010| ) on board the Herschel Space Observatory 
in th e Cosmic Evolution Survey (COSMOS; [Scoville et al. ' 
2007[ ) field. The cataloguej^was constructed by the HerMES 
team, using blind extraction at 250 ^im and a prior-based ex- 
traction for the other two wave bands at 350 /xm and 500 fim. 



Details can be found in Roseboom et al. (20101. Note that 



these sub-millimetre galaxies constitute the lenses here and 
not the sources as in |Wang et al.|p^011| |. 

Unlike 'classical' sub-millimetre sources detected at 
longer wavelength like, e.g., SCUBA-type (Holland et al. 



1999 1 galaxies, the galaxies detected by SPIRE at 250 pm do 



not show a strong negative fc-correction ( [Franceschini et ah] 
1991 Blain & Longair 19931. Hence their redshift distri- 



bution does not extend to as high redshifts, which is benefi- 
cial for the magnification measurement presented here, since 
it allows for easier separation between the sub-millimetre 
lenses and the LEG sources in redshift. 

The beam size of Herschel at 250 ^im is FWHM ^ 18". 
The centroids of the sub-millimetre galaxies in each sky co- 
ordinate are known to Aa = AS = 0.6FWHM/ SNR, where 
SNR is the signal-to- noise ratio of the detection (Ivison et al 



2007[ ). The mean SNR of our lenses is SNR = 5.3, which cor- 
responds to centroid errors of Aa = AS = 2'.'0. Additionally, 
the pointing accuracy of Herschel is limited to ~ 2" as well 
( jPilbratt et al.|[2010[ ). We add these two contributions to 
the centroid error in quadrature and account for this angu- 
lar uncertainty in the modelling (see Sect.|4|. 



2.2 Optical data 

The central part of the COSMOS field (1 deg^) was targeted 
with the MegaCam imager on the Canada France Hawaii 
Telescope (CFHT) as part of the CFHT Legacy Survey. Very 
deep images (5-a limiting magnitudes for point sources of 



26.5-28.0) in the iigriz-bands from the CARS project (Erben 
|et al.|[2009l [Hildebrandt et al.|[2009a| ) were used to select a 
sample of ~ 17 500 Lyman-break galaxies ( [Hildebrandt et ah] 
2009a|b[ ) at redshifts z ~ 3 - 5. These are too distant to be 
resolved in the optical data, but they constitute an ideal 
sample for magnification measurements. 

This set of LBGs has been studied in detail in the lit- 
erature. Their clustering properties are described in [Hilde^ 
brandt et al. (2009a I and a measurement of the LBG lumi- 



nosity function based on this data set is presented in [van 



der Burg et al. (20101. Furthermore, these LBGs were al- 



ready used as background sources for a magnification mea- 
surement in a previous study ( [Hildebrandt et al.|2009b] ) , but 
then using normal galaxies selected by photometric redshifts 
as lenses. 



0.08 



0.06 - 



0.04 



0.02 




^ available at the HeDaM on-line data base: [http: //hedam. oamp. [ 
[fr/HerMES7[ 



Figure 1. Redshift distributions of the sub-millimetre lenses 
(dashed: simulation; dotted: photometric redshifts) and the LBG 
sources (solid) used in the magnification measurement. The sep- 
aration in redshift is not perfect and therefore we model the con- 
tribution from physical clustering (due to redshift overlap) to the 
cross-correlation signal. 



2.3 Redshift distributions 

The HerMES catalogue contains 3402 objects in the common 
area of the optical and sub-millimetre data. For the mass 
measurement we only consider lenses with a sub-millimetre 
fiux density at 250 /xm of 5*250 > 15mJy. This is neces- 
sary because confusion (and other effects) is more severe at 
fainter flux densities in the Herschel data, which would make 
the estimation of a redshift distribution for the lenses very 
problematic. This redshift distribution is needed to interpret 
the signal, since it depends on the lens-source geometry, and 
to estimate a mass. We further apply a colour cut on the sub- 
millimetre lenses to enhance the separation in redshift. We 
require the ratio of the flux densities at 500 [xm and 250 /im 
to be 5500/5*250 < 0.5. This leaves us with 587 lenses. The 
measured and simulated redshift distributions of the lenses 
are presented in [Bethermin et al.[ ( [201l| [2012b[ ). We use the 
simulated one in the following. The model of the redshift 
distribution of the sources is taken from [Hildebrandt et al.l 
( [2009a[ ). These distributions are plotted in Fig. [l] showing 
only marginal overlap. 



3 TECHNIQUE 

The measurement described here is technically similar to 
the previous measurement by [Wang et al.[p011[ ) using sub- 
millimetre galaxies and magnification, but conceptually it 
is very different because here we are studying the effects 
of lensing hy sub-millimetre galaxies rather than the effects 
of lensing by other sources on sub-millimetre galaxies (see 
also [Blake et al.||200^ for a related cross-correlation mea- 
surement). Thus, we can study how their dark matter ha- 
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los deflect light and estimate their average mass through 
this effect, which was not possible when the sub-millimetre 
galaxies were used as sources. 

As detailed in other studies ( Van Waerbeke et al.|20f0 



[Hildebrandt et al.|2011[ ) magnification is a particularly use- 
ful tool for measuring the average masses of objects at high 
redshifts because - unlike shear-based weak lensing methods 
- it does not require the background sources to be resolved. 
A large fraction of the sub-millimetre galaxies studied here 
have redshifts z > 1. With the best optical ground-based 
data there are very few objects at such high redshifts that 
could be resolved and used as background sources for shear 
measurements because of the seeing-limited size of the opti- 
cal point-spread- function (PSF). Even with optical imaging 
data from the Hubble Space Telescope it becomes increas- 
ingly difficult to resolve objects with z > 1.5 due to PSF size. 
Thus, the population of sub- millimetre galaxies studied here 
represents an ideal lens sample for magnification when LBGs 
are used as background sources. 

The WL-magnification effect of the lenses increases the 
fluxes of background objects and shifts their positions on the 
sky. This induces a change in their number density, which 
leads to angular correlations in the positions of lenses and 
sources on the sky, even in the absence of physical cluster- 
ing, i.e., when both populations are well separated in red- 



shift (IScranton et al.||2005l 'Hildebrandt et al.||2009b|). De 



pending on the slope of the magnitude number counts of the 
source sample, a — 2.5 dlog [N{m)]/dm, a positive or a neg- 
ative cross-correlation is expected from WL theory. Bright 
galaxies with typically steep slopes (a — 1 > 0) should show 
positive correlations, and faint galaxies, with shallow slopes 
(a — 1 < 0), should show anti-correlations. 

WL magnification is not the only effect that can cross- 
correlate the sky position of objects that are far apart along 
the line of sight. Extinction by dust in the foreground objects 
can also lead to a characteristic depletion of the number 
density of background objects, and this has to be included 
in the modelling to correctly interpret the cross-correlation 
signal. This behaviour can be described by an effective slope 
aoff as detailed below. 

For this kind of measurement it is particularly impor- 
tant to minimise the overlap in redshift between the lens 
and source populations to suppress a systematic bias due to 
physical cross-correlations, which can be much larger than 
the magnification/extinction signal. 

The cross-correlation functions are estimated using the 
estimator from |Landy fc Szalay| ( [T993[ ) based on pair counts: 



D1D2 - DiR- D2R 
RR 



+ 1, 



(1) 



with D1D2 being the number of sub-millimetre-LBG pairs 
in the angular range [9,9 + 59] normalised by the product 
of their total numbers, DiR being the normalised number 
of pairs between the sub-millimetre/LBG catalogue and a 
random catalogue (with the same surface geometry) in that 
angular range, and RR being the normalised number of pairs 
in the random catalogue in that angular range. By choosing 
a random catalogue that is much larger than the data cat- 
alogues the shot noise introduced by the random catalogue 
is suppressed. 

It has been shown in [Menard et al.| ( |2003[ ) and |Scranton| 
|et al.| ( |2005[ ) that the signal-to-noise ratio of a magnification 



measurement can be boosted when every background galaxy 
is given a weight that corresponds to the a — 1 value at its 
magnitude. The same can be done if there is additional dust 
extinction by instead using an effective weight, a^B — 1. 

For the optimally-weighted correlation function the 
same estimator is used, but D1D2 and D2R are replaced by 



the sums of weights instead of the number of pairs (Hilde- 
[brandt et aL]|2009b| ). The weights are — 1 (see Sect. [5^ 
for a definition) for the LBG in each pair accounting for 
the fact that part of the signal is due to extinction. Such a 
weighting maximises the signal-to-noise ratio of the magnifi- 
cation/extinction measurement ( Menard et al.| [2003 Scran^ 
ton et al.|20 05). The gain is considerable for cases where the 
weight changes appreciably over the whole range of magni- 
tudes of the background sample (typically for very deep data 
as used here; see also Fig. [2]). 



4 THEORETICAL BACKGROUND 

We assume that the halos hosting the sub-millimetre galax- 
ies can be described by a single NFW ( [Navarro et al.|1996[ ) 
halo. This model has two parameters, the mass M200 and the 
concentration c. Detailed dark matter N-body simulations 
show that these two parameters are strongly correlated and 
a relation can be established between the two. Here, we use 
the relation by [Prada et al.| ( [2012[ ) to reduce the number of 
parameters in our mass model. The magnification profile of 
an NFW halo is described in [Wright fc Brainerd| ( |2000| ). 

We model the signal of the optimally-weighted cross- 
correlation function in the following way: 



Uopt 



(61) = {9) + wr (9) + w,, {9) 



(2) 



where w^j describes the contribution from magnification, Wt 
describes the contribution from dust extinction in the lens 
galaxies, and Wcc is the physical clustering part of the signal 
that is due to redshift overlap between lenses and sources 
and should be minimised. 

The magnification signal can be calculated from lensing 
theory: 

"--M (9) = /"""Kh (m) - 1] L (0)"(™)-i - ij iV(m) dm , 

(3) 

being the faintest and brightest magnitudes 



with m 



mm / max 



of the background sample, Qoff(?7i) — 1 being the effective 
weight, n{9) being the magnification profile of an NFW 
halo[^a(m) — 1 = 2.5 dlog [A'^(m)]/dm — 1 being the weight 
in absence of extinction (taken from the LBG luminosity 



van der Burg et al.||2010 l, and N{m) being the 



function of 

normalised number counts of the LBGs. 

The dust absorption A is related to the optical depth 
through 

A = 2.5/ln(10) T = 1.08r . (4) 
We assume that A and the magnification excess Sfi — fi ~ 1 



^ For readability we do not include the redshift dependence here. 
But it should be clear that fJ.{0) has to be calculated by integrating 
over the redshift distributions of the lenses and sources. 
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are related bjj^ 



A — CdSfi . 



(5) 



Here we assume that this is true on average. The redshift 
dependence of the dust extinction and the magnification is 
certainly different. Hence halos at different redshift will con- 
tribute differently to the magnification and extinction sig- 
nals. We neglect this effect in the following but note that 
firm conclusions on the shape of the dust profile can only 
be drawn if such effects are included in the modelling. The 
mass estimate is not directly affected by this simplification. 
Under these assumptions the dust signal becomes: 



™T (9) = / [a^fi (m) - I] 

X [f .08 CdfO"° ''"<™' S^l (9) - f] N{m) dm . (6) 

The contribution from physical clustering to the angular 
correlation function are modelled by 

Wcc{e) = / [aeff (m) - f] bib2WDM{e)N(m) dm , (7) 

where bi/2 are the average bias factors of the lenses and 
sources, respectively, and wdm is the angular correlation 
function of the dark matter fielcQ We conservatively esti- 
mate hi — b2 — 2 in our analysis. 



5 RESULTS 

The angular cross-correlation (unweighted) between the po- 
sitions of all sub-millimetre lenses in the f deg^ area (3402 
objects; no flux cut) and the LBG sources in different mag- 
nitude slices at close separations is shown in Fig. [2] Here we 
use one broad angular bin with f'.'S < 6 < 6'.'0. As expected, 
the positions of the brighter background galaxies are posi- 
tively correlated with the positions of the lenses, while the 
positions of the fainter ones are anti-correlated. However, the 
dashed line, which shows the expected amplitude from mag- 
nification only, without dust extinction, does not fit the data 
points. The solid line represents the expected amplitude of 
the cross-correlation function assuming that on average the 
dust absorption is proportional to the magnification excess 



(note that this is well supported by results from Menard 
|et al.||2010] using SDSS data). We fit for the dust normal- 
isation constant Cd, obtaining a best-fit value of Cd = 0.35. 
The overall normalisation was left as another free parameter 
and both predictions were multiplied by this value. In both 
cases we model N{m) from the LBG luminosity function 



estimates by van der Burg et al. (20f0l 



[Menard et al. (20101 found a value of Cd ~ 0.1 using a 
magnitude-limited, optically-selected, low-z galaxy sample 
from the Sloan Digital Sky Survey. Thus, the sub-millimetre 
galaxies we study here are - not surprisingly - considerably 



^ Note that this does not constrain the dust profile yet because 

Cd can be a function of angular scale. 

* This is calculated with the equation from |Limber| l |l953[ l . We 
used the code by |Hamana et aT] | |2004] l for this. The contribution 
from physical clustering has to be weighed appropriately with the 
product of the redshift distributions of lenses and sources. 
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Figure 2. Amplitude of the angular cross-correlation function be- 
tween 3402 sub- millimetre lenses and LBG sources as a function of 
the LBG i-band magnitude in the angular interval I'.'S < 9 < G'.'O. 
Errors are Poissonian. The dashed line represents the expected 
amplitude if there was only magnification. The solid line repre- 
sents a more realistic model including the effects of extinction by 
dust in the lenses. 



more dusty (although the Cd values cannot be compared di- 
rectly because of the different rest frame wavelengths of the 
filters used for the selection of the lenses). 

Physical clustering signals due to redshift overlap be- 
tween lenses and sources are always positive. Thus, detecting 
a negative cross-correlation, as shown in Fig. [2j already sug- 
gests that we do not have significant redshift overlap here. 

Here, we use all 3402 lenses detected in the sub- 
millimetre data to improve the statistics and be able to 
constrain the scaling between magnification and extinction. 
However, we limit the following analysis to the 587 objects 
with 5250 > 15mJy and S500/S250 < 0.5, for which we can 
estimate a reliable redshift distribution 



(see Sect. [2^. By 
applying the Cd value found from the whole sample of 3402 
objects to the high-confidence sample of 587 objects we im- 
plicitly assume that the redshift distribution, as well as the 
scaling between 5/i and A, does not differ significantly be- 
tween the two samples. With the current data and models 
this cannot be tested directly and has to remain an assump- 
tion here. 

In Fig. [3] the optimally- weighted cross-correlation func- 
tion between these 587 sub-millimetre galaxies and the full 
set of background LBGs is shown, with errors estimated from 
a jack-knife resampling of the background sample. We detect 
a signal at the 7-a level. The best-fit value for Cd (= 0.35) 
that we take from the previous measurement indicates that 
the majority of this signal is due to magnification and not 
extinction. Note that the accurate relative contributions de- 
pend on the absolute value of the magnification itself (and 
hence also on angular scale). At small angular scales mag- 
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Figure 3. Optimally-weighted angular correlation function be- 
tween 587 sub-millimetre galaxies with 250 /im flux density 
5250 > 15 mjy and all background LBGs. The weights are based 
on the results presented in Fig. [5] with only the brightest Lyman- 
break galaxies getting a positive weight and the fainter ones all 
getting a negative weight. The solid line represents the best-fit 
model, consisting of contributions from magnification (dashed 
line), dust extinction (dot-dashed), and a negative contribution 
from physical clustering between lenses and sources (dotted, and 
almost negligible). The top panel shows a logarithmic scale for 
the correlation amplitude whereas the bottom panel has a linear 
scale. Errors are estimated from a jack-knife resampling of the 
background population. 



originating from the measurement of Cd, the assumption that 
dust follows mass, the assumed redshift distributions, or the 
physical clustering model are included here, because these 
uncertainties are very hard to quantify at the current stage, 
but could well be substantial. 

Our best-fit model also suggests that sub-millimetre 
galaxies are very dusty. Assuming that the extinction law 
of the dust in these galaxies is similar to that found in the 
Small Magellanic Cloud ( ^Menard et al.|2010 1, the dust mass 
fraction at the scales we are probing ( ^ 10 kpc) can be esti- 
mated to be roughly Afdust/M2oo « 6 x 10~^. 



6 DISCUSSION 

6.1 Details of the modelling 

The magnification exceeds jj, — 1.5 for the smallest scales 
that we are probing here, taking into account that the dust 
extinction contributes 15% of this signal. Thus, this is not 
strictly weak lensing but rather an intermediate regime be- 
tween strong and weak lensing where a statistical technique 
- similar to what is used in weak lensing - is applied. Note 
that because of this we do not assume the weak lensing limit 
(i.e. Sjj. <C 1) in any part of the modelling. 

Both, the magnification as well as the extinction, lead 
to a net depletion of background sources behind the sub- 
millimetre lenses, hence the negative average weight. Even 
without directly using colour information, we can still dis- 
entangle the two effects because we have a measurement of 
the average intrinsic luminosity function of the LBGs ( |van| 
der Burg et al.|2010 |. Attributing the amplitude of the cor- 
relation function solely to either magnification or extinction 
would lead to inconsistencies as it is shown in Fig. |2] 



nification contributes about 85% to the signal, whereas at 
large angular scales its contribution drops to about 75%. 

We fit the model from Eq. [2] to the data, which repre- 
sents the magnification of an NFW dark matter halo, as well 
as the extinction of a dust halo whose dust is distributed ac- 
cording to the same profile. Additionally we take the small 
contribution from physical cross-correlations w^c into ac- 
count, which is further suppressed because of the weighting 
scheme. Due to the fact that most of the background LBGs 
have negative weights, acff — 1, Wcc actually becomes neg- 
ative in this optimally- weighted casej^We also account for 
errors in the centroid positions of the sub-millimetre lenses 
due to the limited resolution of Herschel. 

The best-fit estimate of the average NFW halo mass 
within r2oo is log^Q [M2oo/Mq] = IS.ntomist&t.) or 
M200 = 1.48l:o 25(stat.) X 10"Mq, directly confirming that 
sub-millimetre galaxies are hosted by very massive halos. 

Note that by assuming that dust follows mass (see 
above) we can exploit the contribution from extinction to the 
signal in order to actually constrain the mass. The quoted 
error contains contributions from shot-noise as well as from 
clustering of the background population. No uncertainties 

^ This also means that not accounting for or underestimating 
uicc would lead to an underestimate of the mass. 



6.2 Relation to Clustering 



Clustering measurements have been used to constrain the 
dark matter halo masses of sub-millimetre galaxies by inter- 
preting their auto-correlation function in the framework of 



the halo model (|Farrah et al. 


2006 


Cooray et al.|2010||Am-| 


blard et al. 2011| Hickox et al. 


2012 


Vieroetal. 20121. While 



we also measure a correlation function here, the origin of the 
signal is very different. The lensing magnification method 
described here differs in some important aspects from such 
clustering measurements. 



We actually try to suppress the physical clustering con- 
tribution to the cross-correlation function as much as possi- 
ble by separating the sub-millimetre lenses and LBG sources 
in redshift. The auto-correlation signal that is used by phys- 
ical clustering measurements represents a systematic nui- 
sance in magnification studies. We model its small contribu- 
tion to the cross-correlation function here and account for it 
in the mass measurement, but it is important to note that 
it is not used to estimate the mass. 

In order to interpret an auto-correlation signal from 
physical clustering, some model for the galaxy bias has to be 
assumed. No such model is needed to interpret our lensing 
correlation function on small scales, under the assumption 
that all sub-millimetre galaxies are central galaxies to their 
halos. 
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Another important difference between physical cluster- 
ing measurements and lensing magnification is the different 
sensitivity to the shape of the redshift distribution. While 
the interpretation of the physical angular cross-correlation 
function depends critically on the width of the redshift dis- 
tribution, lensing is fairly insensitive to this width because 
the lensing efficiency is a slowly changing function of red- 
shift. It is mainly the mean redshift that is important here. 
However, this is only true as long as there is no redshift 
overlap between lenses and sources. Once there are phys- 
ical cross-correlations contributing to the signal, as in the 
measurement presented here, the width becomes important 
again for estimating this contribution. 

6.3 Comparison to abundance matching 

We also roughly estimate the mass of halos hosting the 



sub-millimetre galaxies from abundance matching (Bether- 
mm et al.|[20T2 al) using the following method. We compute 
the infrared luminosity of a source having the mean red- 
shift and the mean flux of the Herschel sample. We convert 
this lu minosity into star formation rate using the |Kennicutt| 
( 1998[) constant (1 x 10'^° A/o/yr/Lo). We finally use the 



abundance matching results at this SFR and this mean red- 
shift, using an interpolation between the two closest redshifts 
where the abundance matching was performed. This yields 
a mass estimate of logjg [M200 / Mq] = 13.2 ± 0.2, in very 
good agreement with our magnification/extinction measure- 
ment rl 



6.4 Notes on the sub-millimetre catalogue 

This sub-millimetre galaxy lens sample represents a fairly 
faint population. Thus, the number of lenses that are lensed 
themselves is also very low. A completely negligible fraction 
of the signal shown in Fig. |3]is due to lensing of both, the 
lenses and the sources, by foreground structures in front of 
both ( [Heavens fc Joachimi|2011| ). 

When we use the full catalogue of 3402 objects for the 
determination of the constant Cd it is very probable that 
there are a large number of spurious objects in this sam- 
ple. Note that this does not matter in that particular case 
because we are only interested in the relative amounts of 
magnification and extinction there. A spurious source should 
not add to either of the two and hence not change the re- 
sult. For the actual mass measurement we use only galaxies 
with 250 /im fiux densities of 5250 > 15mJy as lenses. Given 
the depth of the data this sample should not contain any 
spurious sources. 



6.5 Limitations and possible extensions 

There are several aspects that limit the accuracy of the mea- 
surement presented here. We have to use a redshift distri- 
bution based on a model for the sub-millimetre lenses to 
predict the magnification signal. Individual (photometric) 
redshifts for all sub-millimetre galaxies would certainly help 



^ It should be noted that the estimate from IR abundance match- 
ing is associated with a number of systematic errors that are not 
present in the magnification measurement. 



to overcome any uncertainties in the modelling of their red- 
shift distribution ( [Bethermin et aL]|2012b[ ). 

The amount of physical cross-correlation, described by 
Wcc, depends on the redshift distribution of the LBGs, which 
is taken from simulations (Hildebrandt et al. 2009a I. It 



should be noted that the choice here is actually quite pes- 
simistic. With different plausible choices presented in |Hilde-^ 
brandt et al. (2009a), which are based on modified simula- 
tions, the overlap and hence amplitude of lUcc would decrease 
even further. The mass estimate is virtually unaffected by 
this very small change. We would like to stress that the 
optimal weighting of the correlation function greatly sup- 
presses the physical cross-correlation. In order to boost the 
signal-to-noise ratio of the desired signal components each 
background galaxy is weighted with its expected respon- 
siveness to the combined effects of magnification and ex- 
tinction. These weights are estimated from the background 
galaxies' magnitudes. However, the scaling of the physical 
cross-correlation signal with magnitude of the background 
galaxies is completely different. This suppression is reflected 
in the very low amplitude of Wcc compared to and Wt in 

Fig.m 

We assume that the optical depth of dust extinction in 
the lenses, r, follows the magnification excess, S^, such that 
MK)) ~ CdSfi. This is motivated by the measurement by 



'Menard et al. (20101 with low-z galaxies but has not been 



shown to be true for sub-millimetre galaxies. We can only 
constrain the value of Cd at small scales where the signal-to- 
noise ratio is large enough to split up the LBG background 
sample into several magnitude bins: hence, the choice of the 
broad angular bin with I'.'S < 6 < 6'.'0 in Fig. [2] However, 
our measurement over 1 deg^ is not powerful enough to fully 
constrain the angular dependence of the extinction. Future 
measurements over larger areas could exploit the redden- 
ing effect (Men ard et al.|2010 l to constrain this dependence 
and potentially even the extinction law of the dust in sub- 
millimetre galaxies. 

For consistency, we also calculate models where the dust 
is assumed to be distributed according to an exponential pro- 
file with different scale lengths /ir = 5 — 40kpc instead of 
an NFW-like dust halo. The factor Cd becomes a function of 
angular scale then, and so do the weights for the optimally 
weighted correlation function. The overall amplitude of these 
dust models is fixed by the measured, integrated value of Cd 
at small scales (see Fig. [2|. The best-fit mass for such mod- 
els is smaller by a factor 2-3 than the best-fit mass for the 
model with an NFW-like dust halo - the larger the scale 
length the smaller the discrepancy. However, the reduced 
of these model fits to the data, x^/dof = 2-5, is consider- 
ably larger than the fiducial x^/dof — 1.5 of the model with 
an NFW-like dust component - again with the largest dust 
scale length (40kpc) yielding the smallest x^/dof. This in- 
dicates that the dust is indeed very widely distributed, and 
our measurements favour the fiducial model where "dust fol- 
lows mass". It should be stressed that we do not interpret 
the results in such a way that these sub-millimetre galaxies 
have a smooth NFW-like dust halo or exponential dust disks 
with very large scale-lengths of /ia>40kpc. It is, however, 
possible that some additional, widely extended dust compo- 
nent - similar to what is found in Menard et al. (20101 - is 



responsible for the effect seen here. 

Another consequence of the limited signal-to-noise ratio 
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is that we cannot fit a multi-parameter model to the data 
points in Fig[3] For this reason we concentrate on small an- 
gular scales, where a single halo can be assumed to domi- 
nate (1-halo term). The contributions from other halos (2- 
halo term) only become important at larger angular scales 
and can be neglected here for simplicity. Further assuming a 
mass-concentration relation ( jPrada et al.|2012[ ) for the NFW 
halo leaves us with just one parameter to fit, the average 
mass M200- We check the influence of our choice of a par- 
ticular mass-concentration relation by also implementing the 
relation by Munoz-Cuartas et al. ( 2011[ ) and find a negligible 
difference in the best-fit mass of Alog^g [^^200/^^0] « 0.03. 
Again, the interpretation of a measurement from a more 
powerful, larger-area survey could easily be extended to con- 
strain the halo-occupation distributions, satellite fractions, 
and concentration parameters directly. 

We also assume that all sub-millimetre galaxies are cen- 
tral galaxies to their halos. There are some indications in the 
literature from clustering studies that a fair fraction of the 
galaxies in such a sub-millimetre sample are in fact satellites 
(~ 25% for our lens sample; see Cooray et al.|20T0 l. In the 
regime probed here the magnification excess, 5^, is fairly 
linear in the mass. Thus, assuming a worst-case scenario, 
where no mass would be associated with 25% of our lenses, 
this would lead to an underestimate of the halo mass of the 
remaining 75% central galaxies of AlogjQ [M2oo/-^^©l ~ O-l- 
This is similar to the statistical error of our measurement 
so that we decided not to correct for this directly here, es- 
pecially because the real effect is certainly smaller than this 
worst-case scenario. 



7 SUMMARY AND OUTLOOK 

In this paper we show how to use the magnification ef- 
fect of WL to measure the average mass of sub-millimetre 
galaxies. Using sub-millimetre galaxies as the lenses and 
Lyman-break galaxies as the sources we find a mass of 
log^o [M2OO/M0] = 13.17to og(stat.) for the halos host- 
ing the sub-millimetre galaxies. The presence of significant 
amounts of extinction by dust in the lenses complicates this 
measurement. However, accounting for the dust allows us to 
simultaneously constrain the dust mass fraction of the lenses 
to Mdust/M2oo ~ 6 X 10"^ 

With deep, large-area, imaging surveys (DES, LSST, 
Euclid) on the horizon, WL-magnification methods will gain 
additional importance. The higher the redshifts of the ob- 
jects under study, the fewer the techniques that can pro- 
vide a reliable mass measurement. As we have shown in 
this paper, magnification can provide such mass estimates 
for high redshift objects. Follow-up observations of the opti- 
cal surveys mentioned above with sub-millimetre telescopes 
will provide very large samples of high-z star-forming galax- 
ies. Refined magnification measurements of the kind pre- 
sented here, alongside with clustering measurements and 
other techniques, will yield unprecedented insights into the 
physical processes of star-formation in the first half of the 
Universe. 

Future measurements with better statistics will enable 
the lensing and extinction effects to be separated. Addi- 
tional redshift information for the sub-millimetre galaxies 
will lead to more accurate mass estimates. By splitting the 



sub-millimetre galaxy sample in brightness it will be pos- 
sible to directly study the relationship between mass and 
star-formation rate. 

This study shows that extinction can not be neglected 
in magnification studies. The studies by [Hildebrandt et al.| 
(2009b I and Morrison et al. (20121 show excessive anti- 



correlations when very faint sources are used. Also the am- 
plitudes of the angular correlation function presented there 
turn over from positive to negative at brighter source magni- 
tudes than expected from magnification alone. This hints at 
extinction playing a role. The framework outlined here can 
explain this effect, and we strongly suggest future magnifica- 
tion studies to account for extinction to avoid a systematic 
bias. 



ACKNOWLEDGEMENTS 

Herschel is an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator con- 
sortia and with important participation from NASA. 
SPIRE has been developed by a consortium of institutes led 
by Cardiff University (UK) and including: Univ. Lethbridge 
(Canada); NAOC (China); CEA, LAM (France); IFSI, Univ. 
Padua (Italy); lAC (Spain); Stockholm Observatory (Swe- 
den); Imperial College London, RAL, UCL-MSSL, UKATC, 
Univ. Sussex (UK); and Caltech/JPL, IPAC, Univ. Colorado 
(USA). This development has been supported by national 
funding agencies: CSA (Canada); NAOC (China); CEA, 
CNES, CNRS (France); ASI (Italy); MCINN (Spain); SNSB 
(Sweden); STFC, UKSA (UK); and NASA (USA), 
This research has made use of data from the HerMES 
project (http : //hermes . Sussex . ac.uk/). HerMES is a Her- 
schel Key Programme utilising Guaranteed Time from the 
SPIRE instrument team, ESAC scientists and a mission sci- 
entist. 

The HerMES data were accessed through the HeDaM data- 
base (http://hedaiii.oamp.fr) operated by CeSAM and 
hosted by the Laboratoire d'Astrophysique de Marseille. 
The optical data are based on observations obtained 
with MegaPrime/MegaCam, a joint project of CFHT and 
CEA/DAPNIA, at the Canada-France-Hawaii Telescope 
(CFHT) which is operated by the National Research Coun- 
cil (NRC) of Canada, the Institut National des Sciences de 
rUnivers of the Centre National de la Recherche Scientifique 
(CNRS) of France, and the University of Hawaii. This work 
is based in part on data products produced at TERAPIX 
and the Canadian Astronomy Data Centre, as part of the 
Canada-France-Hawaii Telescope Legacy Survey, a collabo- 
rative project of NRC and CNRS. 

HH is supported by the Marie Curie lOF 252760, a CITA 
National Fellowship, and the DEC grant Hi 1495/2-1. LVW 
is supported by NSERC and ClfAR. DS acknowledges sup- 
port by NSERC and CSA. TE is supported by the Deutsche 
Forschungsgemeinschaft through project ER 327/3-1 and 
the Transregional Collaborative Research Centre TR 33 - 
aAlThe Dark UniverseaAl. RVDB acknowledges support 
from the Netherlands Organisation for Scientic Research 
grant number 639.042.814. 
Author Contributions: 

H.H. led the analysis and wrote the draft version of the pa- 
per. L.v.W. and D.S. contributed significantly in the development 



© 2011 RAS, MNRAS 000,[T}j9] 



Inferring the mass of sub-millimetre galaxies through magnification 9 



of this project through ideas and discussions on a daily basis. M.B. 
provided the redshift distribution of the sub-milUmetre galaxies. 
T.E. led the optical data reduction and calibration. R.F.J.v.d.B. 
provided the luminosity function estimate of the background 
galaxies. All other co-authors contributed extensively and equally 
by their varied contributions to the SPIRE instrument, the Her- 
schel mission, analysis of SPIRE and HcrMES data, planning of 
HerMES observations and scientific support of HcrMES, and by 
commenting on this manuscript as part of an internal review pro- 
cess. 



REFERENCES 

Amblard A. et al., 2011, Nature, 470, 510 

Bartelmann M., Schneider P., 2001, Phys. Rep., 340, 291 

Bethermin M., Dole H., Lagache G., Le Borgne D., Penin 

A., 2011, A&A, 529, A4 
Bethermin M., Dore O., Lagache G., 2012a, A&A, 537, L5 
Bethermin M. et al., 2012b, arXiv:1203.1925 
Blain A. W., Longair M. S., 1993, MNRAS, 264, 509 
Blake C., Pope A., Scott D., Mobasher B., 2006, MNRAS, 

368, 732 

Cooray A. et al., 2010, A&A, 518, L22 
Eggen O. J., Lynden-BeU D., Sandage A. R., 1962, ApJ, 
136, 748 

Erben T. et al., 2009, A&A, 493, 1197 

Faber S. M., Jackson R. E., 1976, ApJ, 204, 668 

Farrah D. et al., 2006, ApJ, 641, L17 

Ferrarese L., Merritt D., 2000, ApJ, 539, L9 

Franceschini A., Toffolatti L., Mazzei P., Danese L., de 

Zotti G., 1991, A&AS, 89, 285 
Griffin M. J. et al., 2010, A&A, 518, L3 
Hamana T., Ouchi M., Shimasaku K., Kayo I., Suto Y., 

2004, MNRAS, 347, 813 
Heavens A. F., Joachimi B., 2011, MNRAS, 415, 1681 
HerMES Collaboration et al., 2012, ArXiv e-prints 
Hickox R. C. et al., 2012, MNRAS, 421, 284 
Hildebrandt H. et al., 2011, ApJ, 733, L30-I- 
Hildebrandt H., Pielorz J., Erben T., van Waerbeke L., 

Simon P., Capak P., 2009a, A&A, 498, 725 
Hildebrandt H., van Waerbeke L., Erben T., 2009b, A&A, 

507, 683 

Holland W. S. et al., 1999, MNRAS, 303, 659 
Ivison R. J. et al., 2007, MNRAS, 380, 199 
Kennicutt, Jr. R. C., 1998, ApJ, 498, 541 
Komatsu E. et al., 2011, ApJS, 192, 18 
Landy S. D., Szalay A. S., 1993, ApJ, 412, 64 
Leauthaud A. et al., 2012, ApJ, 744, 159 
Limber D. N., 1953, ApJ, 117, 134 

Madau P., Ferguson H. C., Dickinson M. E., Giavalisco M., 

Steidel C. C., Fruchter A., 1996, MNRAS, 283, 1388 
Menard B., Hamana T., Bartelmann M., Yoshida N., 2003, 

A&A, 403, 817 
Menard B., Scranton R., Fukugita M., Richards G., 2010, 

MNRAS, 405, 1025 
Morrison C. B., Scranton R., Menard B., Schmidt S. J., 

Tyson J. A., Ryan R., Choi A., Wittman D. M., 2012, 

MNRAS, 426, 2489 
Muiioz-Cuartas J. C, Maccio A. V., Gottlober S., Dutton 

A. A., 2011, MNRAS, 411, 584 
Navarro J. F., Frenk C. S., White S. D. M., 1996, ApJ, 462, 

563 



Oliver S. J. et al., 2010, A&A, 518, L21 

Pilbratt G. L. et al., 2010, A&A, 518, LI 

Prada F., Klypin A. A., Cuesta A. J., Betancort-Rijo J. E., 

Primack J., 2012, MNRAS, 423, 3018 
Roseboom I. G. et al., 2010, MNRAS, 409, 48 
ScoviUe N. et al., 2007, ApJS, 172, 1 
Scranton R. et al., 2005, ApJ, 633, 589 
Steidel C. C, Giavalisco M., Pettini M., Dickinson M., 

Adelberger K. L., 1996, ApJ, 462, L17 
TuUy R. B., Fisher J. R., 1977, A&A, 54, 661 
van der Burg R. F. J., Hildebrandt H., Erben T., 2010, 

A&A, 523, A74 
Van Waerbeke L., Hildebrandt H., Ford J., Milkeraitis M., 

2010, ApJ, 723, L13 
Viero M. P. et al., 2012, ArXiv e-prints 
Wang L. et al, 2011, MNRAS, 414, 596 
Wright C. O., Brainerd T. G., 2000, ApJ, 534, 34 



© 2011 RAS, MNRAS 000,[T]{9] 



